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Abstract 



Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized 
with MPI, we investigated long-term particle acceleration associated with a relativis- 
tic electron-positron jet propagating in an unmagnetized ambient electron-positron 
plasma. The simulations were performed using a much longer simulation system 
than our previous simulations in order to investigate the full nonlinear stage of the 
Weibel instability and its particle acceleration mechanism. Cold jet electrons are 
thermalized and ambient electrons are accelerated in the resulting shocks. Accelera- 
tion of ambient electrons leads to a maximum ambient electron density three times 
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larger than the original value. Behind the bow shock in the jet shock strong elec- 
tromagnetic fields are generated. These fields may lead to time dependent afterglow 
emission. We calculated radiation from electrons propagating in a uniform parallel 
magnetic field to verify the technique. We also used the new technique to calculate 
emission from electrons based on simulations with a small system. We obtained 
spectra which are consistent with those generated from electrons propagating in 
turbulent magnetic fields with red noise. This turbulent magnetic field is similar 
to the magnetic field generated at an early nonlinear stage of the Weibel instabil- 
ity. A fully developed shock within a larger system generates a jitter /synchrotron 
spectrum. 

Key words: acceleration of particles, galaxies, jets, gamma rays bursts, magnetic 
fields, plasmas, shock waves, radiation 



1 RPIC Simulations 



Particle-in-cell (PIC) simulations can shed light on the physical mechanism 
of particle acceleration that occurs in the complicated dynamics within rela- 
tivistic shocks. Recent PIC simulations of relativistic electron-ion and electron- 
positron jets injected into an ambient plasma show that acceleration occurs 
within the downstream jet . In general, these simulations 

have confirmed that relativistic jets excite the Weibel instability, which gen- 



erates current filaments and associated magnetic fields (jlOl ). and accelerates 
electrons fiU § 8 3 B 0) ■ 

Therefore, the investigation of radiation resulting from accelerated particles 
(mainly electrons and positrons) in turbulent magnetic fields is essential for 
understanding radiation mechanisms and their observable spectral properties. 
In this report we present a new numerical method to obtain spectra from 
particles self-consistently traced in our PIC simulations. 

Pair Jets Injected into Unmagnetized Pair Plasmas using a Large 

System 

We have performed simulations using a system with (L x , L y , L z ) = (4005A, 
131A, 131A) (A = 1: grid size) and a total of ~ 1 billion particles (12 
particles/cell/species for the ambient plasma) in the active grid zones (llll ). 
In the simulations the electron skin depth, A ce = c/u pe = 10. OA, where 
Lu pc = (Ane 2 n e /m e ) l l 2 is the electron plasma frequency and the electron Debye 
length A e is half of the grid size. Here the computational domain is six times 
longer than in our previous simulations (j3;0). The electron number density of 
the jet is 0.676n e , where n e is the ambient electron density and 7 = 15. The 
electron/positron thermal velocity of the jet is v? th = 0.014c, where c = 1 is 



2 



60 




500 1000 1500 2000 2500 3000 3500 



Fig. 1. The averaged values of electron density (a) and field energy (b) along the x 
at t = 3750WPC 1 . Fig. la shows jet electrons (red), ambient electrons (blue), and the 
total electron density (black). Fig. lb shows electric field energy (red) and magnetic 
field energy (blue) divided by the total kinetic energy. 
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Fig. 2. The case with a strong magnetic field (-E> x = 3.7) and larger perpendicular 
velocity (v±i = 0.1c, v±2 = 0.12c). The paths of two electrons moving helically along 
the x— direction in a homogenous magnetic field shown in the x — y-plane (a). The 
two electrons radiate a time dependent electric field. An observer situated at great 
distance along the n-vector sees the retarded electric field from the moving electrons 
at the rest frame (b). The observed power spectrum at different viewing angles from 
the two electrons (c). Frequency is in units of u>~ e . 

the speed of light. 

Figure 1 shows the averaged (in the y — z plane) electron density and elec- 
tromagnetic field energy along the jet at 3750co>~ e 1 . The resulting profiles of 
jet (red), ambient (blue), and total (black) electron density are shown in Fig. 
la. The ambient electrons are accelerated by the jet electrons and pile up 
towards the front part of jet. At the earlier time the ambient plasma density 
increases linearly behind the jet front. At the later time the ambient plasma 
shows a rapid increase to a plateau behind the jet front, with additional in- 
crease to a higher plateau farther behind the jet front. The jet density remains 
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approximately constant except near the jet front. 
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Fig. 3. Two-dimensional images in the x — z plane at y/A = 65 for t = 450u;p e 1 . 
The colors indicate the x-component of current density generated by the Weibel 
instability, with the x- and z-components of magnetic field represented by arrows 
(a). Phase space distributions as a function of x/A — *yv x plotted for the jet (red) 
and ambient (blue) electrons at the same time. 

The Weibel instability remains excited by continuously injected jet particles 
and the electromagnetic fields are maintained at a high level, about four times 
that seen in a previous, much shorter grid simulation system (L x = 640 A). At 
the earlier simulation time a large electromagnetic structure is generated and 
accelerates the ambient plasma. As shown in Fig. lb, at the later simulation 
time the strong magnetic field extends up to x/A = 2, 000. These strong fields 
become very small beyond x/A — 2000 in the shocked ambient region (0; 0). 

The acceleration of ambient electrons becomes visible when jet electrons 
pass about x/A = 500. The maximum density of accelerated ambient electrons 
is attained at t — 1750^"^. The maximum density gradually reaches a plateau 
as seen in Fig. la. The maximum electromagnetic field energy is located at 
x/A = 1, 700 as shown in Fig. lb. 



1.1 New Numerical Method for Calculating Synchrotron Emission 



Let a particle be at position r (t) at time t (1121 ; Il3l ; Il4h . At the same 



time, we observe the electric field from the particle from position r. However, 
because of the finite velocity of light, we observe the particle at an earlier 
position r (t ) where it was at the retarded time t — t — 8t = t — R(t )/c. 
Here R(t ) = |r — r (t )| is the distance from the charge (at the retarded time 
t ) to the observer. 
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Fig. 4. Two-dimensional images in the x — z plane at y/A = 65 for t = 525cjp e . 
The colors indicate the x-component of current density generated by the Weibel 
instability, with the x- and z-components of magnetic field represented by arrows 
(a). Phase space distributions as a function of x/A — *yv x plotted for the jet (red) 
and ambient (blue) electrons at the same time. 

After some calculation and simplifying assumptions the total energy W 
radiated per unit solid angle per unit frequency from a charged particle mo ving 
with instantaneous velocity (3 under acceleration can be expressed as (15ul6l) 



d 2 W 



/j cq 



dVtduj 167T 3 



n x [(n - 0) x $] 



,iu(t -nr (t )/c) ( j£ 



n 



Here, n = R(t')/|R(t')| is a unit vector that points from the particle's 
retarded position towards the observer. 

The observer's viewing angle is set by the choice of n (n^+n^+nl = 1). The 
choice of unit vector n along the direction of propagation of the jet (hereafter 
taken to be the x-axis) corresponds to head-on emission. For any other choice 
of n (e.g., # 7 = I/7), off-axis emission is seen by the observer. 

In order to calculate radiation from relativistic jets propagating along the x 
direction (11 21 ) we consider a test case which includes a parallel magnetic field 
(£> x ), and jet velocity of fji^ = 0.99c. Two electrons are injected with different 
perpendicular velocities (v±i = 0.1c, v± 2 = 0.12c). A maximum Lorenz factor 
of 7max = {(1 — (vf 2 + ^Vc 2 } -1 ^ 2 = 13.48 is calculated with the larger 



perpendicular velocity. 

Figure 2 shows electron trajectories in the x 



y plane (red: v± 2 = 0.12c, 
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blue: v±\ = 0.1c) (a: left panel), the radiation (retarded) electric field (b: mid- 
dle panel), and spectra (right panel) for the case B x = 3.70. The two electrons 
are propagating left to right with gyration in the y — z plane (not shown). 
The gyroradius is about 0.44A for the electron with the larger perpendicu- 
lar velocity. The seven curves show the power spectrum at viewing angles of 
0° (red), 10° (orange), 20° (yellow), 30° (moss green), 45° (green), 70° (light 
blue), and 90° (blue). The higher frequencies become stronger at the 10° view- 
ing angle. The critical angle for off-axis radiation # 7 = 180°/ (7T7 max ) for this 
case is 4.25°. As shown in this panel, the spectrum at a larger viewing angle 
(> 20°) has smaller amplitude. 



Since the jet plasma has a large velocity x-component in the simulation 
frame, the radiation from the particles (electrons and positrons) is strongly 



beamed along the x-axis (jitter radiation) (1171 ; Il8l ) 



Equations 6.30a and 6.30b show that the radiation with the viewing ang le 
a = disappears (see Fig. 6.5 in the textbook of Rybicki and Lightman (Il5l)). 
However, based on two other textbooks, radiation at the viewing angle 0° 
should not vanish (Il9t |20j). This aspect is shown in Fig. 2c, and at the higher 
frequency the amplitude at the viewing angle 10° is stronger than that with 
viewing angle 0°. 



1.2 The Standard Synchrotron Radiation Model 



A synchrotron shock model is widely adopted to describe the radiation 
mechanism in the external shock thought to be responsible for observed broad- 
band GRB afterglows Q B H 0) ■ Associated with this model are three 
major assumptions that are adopted in almost all current GRB afterglow mod- 
els. Firstly, electrons are assumed to be "Fermi" accelerated at the relativistic 
shocks and to have a power-law distribution with a power-law index p upon 
acceleration, i.e. N(E e )dE c oc E~ p dE c . This is consistent with recent PIC 
simulations of the shock formation and particle acceleration (Q) and also some 
Monte Carlo models 800), but see (@;[3o|). Secondly, a fraction e P (gen- 
erally taken to be < 1) of the total electrons associated with ISM baryons are 
accelerated, and the total electron energy is a fraction e e of the total internal 
energy in the shocked region. Thirdly, the strength of the magnetic fields in 
the shocked region is unknown, but its energy density (_B 2 /87r) is assumed to 
be a fraction of the internal energy. These assumed "micro-physics" pa- 



rameters, p, e e and 6b, whose values are obtained from spectral fits (13ll ; l32l ) 



reflect a lack of knowledge of the underlying microphysics (1331 ) . 



The typical observed emission frequency from an electron with (comoving) 
energy 7 e m e c 2 in a frame with a bulk Lorentz factor r is v = (eB/27rm e c). 
Three critical frequencies are defined by three characteristic electron energies. 
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These are is m (the injection frequency), v c (the cooling frequency), and um (the 
maximum synchrotron frequency). In our simulations of GRB afterglows, there 
is one additional relevant frequency, i/ a , due to synchrotron self-absorption at 



lower frequencies (34;S; 24; 25h . 



The general agreement between the blast wave dynamics and the direct 
measurements of the fireball size argue for the validity of this model's dynamics 
( 24 ; 1251) . The shock is most likely collisionless, i.e. mediated by plasma instabil- 
ities (1331 ). The electromagnetic instabilities mediating the afterglow shock are 
expected to generate magnetic fields. Afterglow radiation was therefore pre- 
dicted to result from synchrotron emission of shock accelerated electrons (1361 ). 
The observed spectrum of afterglow radiation is indeed remarkably consistent 
with synchrotron emission of electrons accelerated to a power-law distribu- 
tion, providing support for the standard afterglow model based on synchrotron 



emission of shock accelerated electrons (J39|; |40|; |22j; |2l|; |37J; |38|; [2J; |25| ) 



In order to determine the luminosity and spectrum of synchrotron radia- 
tion, the strength of the magnetic field (ee) and the energy distribution of the 
electrons (p) must be determined. Due to the lack of a first principles the- 
ory of collisionless shocks, a purely phenomenological approach to the model 
of afterglow radiation was ascribed without investigating in detail the pro- 
cesses responsible for particle acceleration and magnetic field generation (1331 ). 
Rather, one simply assumes that a fraction e-Q of the post-shock thermal en- 
ergy density is carried by the magnetic field, that a fraction e e is carried by 
electrons, and that the energy distribution of the electrons is a power-law, 
d log n e /d log e = p (above some minimum energy Eq which is determined by e e 
and p), 6b, e e and p are treated as free parameters, determined by observations. 
It is important to clarify here that the constraints implied on these parameters 
by the observations are independent of any assumptions regarding the nature 
of the afterglow shock and the processes responsible for particle acceleration 
or magnetic field generation. Any model should satisfy these observational 
constraints. 



The properties of synchrotron (or "jitter") emission from relativistic shocks 
will be determined by the magnetic field strength and structure and the elec- 
tron energy distribution behind the shock. The characteristics of jitter radia- 
tion may be important to understanding the complex time evolution and/or 
spectral structure in gamma-ray bursts (J4lj). For example, jitter radiation 
has been proposed as a means to explain GRB spectra below the peak fre- 
quency that are harder than the "line of death" spectral index associated 
with synchrotron emission ( 17c Il8l ). i.e., the observed spectral power scales as 
F u oc z/ 2//3 , whereas synchrotron spectra are F u oc v 1 ^ or softer (fl8l ). Thus, it 
is essential to calculate radiation production by tracing electrons (positrons) 
in self-consistently treated small-scale electromagnetic fields. 
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1.3 Calculating Synchrotron and Jitter Emission from Electron Trajectories 
in Self- consistently Generated Magnetic Fields 

In order to obtain the spectrum of synchrotron (jitter) emission, we consider 
an ensemble of electrons selected in the region where the Weibel instability 
has fully grown and electrons are accelerated in the generated magnetic fields. 
In order to validate our numerical method we performed simulations using 
a small system with (L x , L y , L z ) = (645A, 131A, 131A) (A = 1: grid size) 
and a total of ~ 0.5 billion particles (12 particles/cell/species for the ambient 
plasma) in the active grid zones (0). First we performed simulations without 
calculating radiation up to t = 450^"^. The jet front is located around about 
x/A = 480. We selected 12,150 electrons for each jet and ambient electrons 
randomly. Recently, a similar calculation has been carried out for the radiation 
from accelerated electrons in laser-wakefield acceleration (j42l ) and in shocks 

a. 

Figure 3 shows (a) the current filaments generated by the Weibel instability 
and (b) the phase space of x/A — jV x for jet electrons (red) and ambient 
electrons (blue) at t = 450a;p e 1 . 

Iog, f((u) 
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Fig. 5. Spectra obtained from jet and ambient electrons for the two viewing angles. 
Spectra with jet electrons are shown in red (0°) and orange (5°). Spectra from 
ambient electrons show the lowest levels by blue (0°) and light blue (5°). 

Figure 4 shows (a) the x-component of current density generated by the 
Weibel instability and (b) the phase space of jet electrons and ambient elec- 
trons after t s = 75u~^ (at t = 525co~ e 1 ). 

We calculated the emission from 12,150 electrons during the sampling time 
t B = t 2 — ti = 75o»p e 1 with Nyquist frequency = 1/2 At = 200co> pe where 
At = O.OOSu;" 1 is the simulation time step and the frequency resolution Auj = 
l/t s = 0.0133^. 

The spectra shown in Fig. 5 are obtained for emission from jet electrons and 
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Fig. 6. The case with a weak magnetic field (B x = 0.37) and small perpendicular 
velocity (v±i = 0.01c, v±2 = 0.012c). The paths of two electrons moving helically 
along the x— direction in a homogenous magnetic field shown in the x — y-plane (a). 
The two electrons radiate a time dependent electric field. An observer situated at 
great distance along the n-vector sees the retarded electric field from the moving 
electrons (b). The observed power spectrum at different viewing angles from the 
two electrons (c) . Frequency is in units of uj^ . 

ambient electrons separately. In this case the spectra are calculated for head-on 
radiation (0°) and 5°. The radiation from jet electrons show Bremsstrahlung- 



like spectra as a red line (0°) and orange line (5°) (1131 ). The spectra with jet 
electrons are different from the spectra shown in Fig. 2c. Since the magnetic 
fields generated by the Weibel instability are rather weak and the jet electrons 
are not much accelerated, the trajectories of jet electrons are almost straight 
with only a slight bent. 

We compare these spectra with our known spectra obtained from two (jet) 
electrons, the case with a parallel magnetic field (B x = 0.37), and jet velocity of 
v ji,2 — 0.99c. Two electrons are injected with different perpendicular velocities 
(v±i = 0.01c, v± 2 = 0.012c). A maximum Lorenz factor of 7 max = {(1 — 
(v? 2 + Vj_ 2 ) I 'c 2 } -1 / 2 = 7.114 accompanies the larger perpendicular velocity. 
The critical angle for off-axis radiation 9 1 = 180°/ (7T7 max ) for this case is 
8.05°. 



Comparing the spectra with Figs. 5 and 6c we find similarities. The lower 
frequencies have flat spectra and the higher frequencies decrease monotoni- 
cally. The slope in Fig. 5 is less steep than that in Fig. 6c. This is due to the 
fact that the spread of Lorenz factors of jet electrons is larger and the average 
Lorenz factor is larger as well. Furthermore, even the magnetic field strength is 
not so large, however the slope of the spectra seems to be consistent with Fig. 
7.16 (left) with the turbulent magnetic field with the red noise (/x = —3) in 
Hededal's Ph. D. thesis (Il3l ). We obtained several different parameters with jet 
electrons and ambient magnetic field. However, the strength of the magnetic 
fields generated by the Weibel instability is small, therefore the spectra for 
these cases are very similar to Fig. 5. As shown in Fig. 7.12 in Hededal's Ph. 
D. thesis (flil ). the trajectories of jet electrons have to be chaotic to produce 
a jitter- like spectrum Fig. 7.22 (1131 ). 
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As shown in Fig. lb, the magnetic field energy in the region x/A< 500 is 
small (cb < 0.07), therefore, as expected, the spectra look like those emitted 
from electrons propagating in a turbulent magnetic field with red noise (see 
also Figs. 3a and 4a). 



1.4 Discussions 



Emission obtained with the method described above is obtained self-consistently, 
and automatically accounts for magnetic field structures on small scales re- 
sponsible for jitter emission. By performing such calculations for simulations 
with different parameters, we can investigate and compare the different regimes 



of jitter- and synchrotron-type emission (1171 ; Il8l ). The feasibility of this ap 



proach has already been demonstrated (Il3l ; Il4j ). and its implementation is 



straightforward. Thus, we should be able to address the low frequency GRB 



spectral index violation of the synchrotron spectrum line of death (Il8l ). 



Medvedev and Spitkovsky recently showed that electrons may cool effi- 
ciently at or near the shock jump and are capable of emitting a large fraction 



of the shock energy (1441 ) . Such shocks are well- resolved in existing PIC sim- 
ulations; therefore, the microscopic structure can be studied in detail. Since 
most of the emission in such shocks would originate from the vicinity of the 
shock, the spectral power of the emitted radiation can be directly obtained 
from finite-length simulations and compared with observational data. 

As shown in Fig. 1, behind the trailing shock the electrons are accelerated 
and strong magnetic fields are generated. Therefore, this region seems to pro- 
duce the emission that is observed by satellites. We will calculate more spectra 
based on our RPIC simulations and compare in detail with Fermi data. 
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